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ABSTRACT 

We investigate the luminosity functions (LFs) and projected number density profiles 
of galactic satellites around isolated primaries of different luminosity. We measure 
these quantities for model satellites placed into the Millennium and Millennium II 
dark matter simulations by the GALFORM semi-analytic galaxy formation model 
for different bins of primary galaxy magnitude and we investigate their dependence 
on satellite luminosity. We compare our model predictions to the data of Guo et 
al. from the Sloan Digital Sky Survey Data Release 8 (SDSS DR8). First, we use a 
mock light-cone catalogue to verify that the method we used to count satellites in 
the SDSS DR8 is unbiased. We find that the radial distributions of model satellites 
are similar to those around comparable primary galaxies in the SDSS DR8, with 
only slight differences at low luminosities and small projected radii. However, when 
splitting the satellites by colour, the model and SDSS satellite systems no longer 
resemble one another, with many red model satellites, in contrast to the dominant 
blue fraction at similar luminosity in the SDSS. The few model blue satellites are 
also significantly less centrally concentrated in the halo of their stacked primary than 
their SDSS counterparts. The implications of this result for the GALFORM model 
are discussed. 

Key words: galaxies: dwarf, galaxies: structure, (galaxies:) Local Group, Galaxies: 
fundamental parameters, galaxies: statistics 



INTRODUCTION 



While the standard A Cold Dark Matter (ACDM) model has 
been shown to be in good agreement with observations of 
large scale structure, the verification of this model at small, 
galactic scales remains less certain. One reason for this is 
the increased importance of astrophysical processes relative 
to gravity in this strongly non-linear regime. The study of 
the properties and distribution of galactic satellite galaxies 
provides an opportunity to test ACDM on small scales while 
also constraining different aspects of galaxy formation mod- 
els related to the rates at which satellites form stars, become 
disrupted and merge with the central galaxy. 

The Local Group satellite system within the haloes of 
the Milky Way (MW) and M31 is often the focus of studies 
of galactic satellites because it is here where the lowest lumi- 
nosity satellites can be detected. While this system has been 
used in attempts t o constrain the cosmological model (e 
Bullock et al.ll2000l: iBenson et al]|2002i ; iKlypin et al.ll200 
Lovell et al.ll2012l ; IWang et al.ll2012l ). it is not clear that it is 
typical of the population as a whole. Using the technique of 
abundance matching in the Millennium II simulation (MS- 
II), a large- volume, high-res olution dark matter simulation, 
iBovlan-Kolchin et al.l feoid ) concluded that there should be 
significant scatter in the properties of satellite systems from 



one prima r y to a nother. This conclusion was supported by 
ICuo et all (|201ll ) who applied a semi-analytic galaxy for- 
mation model to the same simulation. Thus, a large, statis- 
tically representative sample of primary galaxies is clearly 
needed to test cosmological and galaxy formation models. 
Such an approach also avoids the difficulty of having to de- 
fine quite what is meant by a Local Group. 

The construction of large galaxy redshift surveys, such 
as the Two Degree Field Gala xy Redshift Surve y (2dFGRS, 
IColless et ai1l200ll ). the SDSS dYork et alfeOOOT), the G alaxy 
and Mass Assembly (GAMA, iDriver et al.ll2009l . 1201 ll ) Sur- 
vey and the Deep Extra galactic Evolutionary Probe 2 
(DEEP2, iDavis et al1l2003l ) Survey, has led to the accumu- 
lation of external galaxy samples covering large volumes. 
Many studies of the luminosity function (LF), spatial distri- 
bution and kinematics of bright satellites ha ve been carried 
out using these and even earlier data sets (e.g. Zaritskv et al.l 
1991 Il997l: Ivan den Bosch et al.1 12004 IConrov et al.l 120051: 



van den Bosch et al. 1120051: IChen et alll2006l: [Koposov et al.l 
20091 : iBusha et al.l I2010J : IPrescott et al.l 1201 if ). The inchT 
sion of a statistical background subtraction in the satel- 
lite system estimators has allowed fainter satellites to be 
studied using deeper photometric galaxy catalogues (e.g 
Lorrimer et al.lll994l;lLiu et al.ll2008l. l201ll: [Lares et alj|2011 



Guo et al.ll201ll : lNierenberg et al.ll201ll : lTal fc van Dokkuml 
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the regime over which the satellite distributions have been 
quantified, a more stringent test of the models can be per- 
formed. 

While some studies have attempted to make model 
predictions u sing cosmologi c al h yd rodynamic simu- 
lations (e.g Libeskind et al.l 2007 ; Okamoto fe Frenkl 



20091; lOkamoto et all l20ld: IWadepuhl fe Springell l201ll: 



Parry et all l2012h . these efforts are limited to very few 



primary galaxies because of the high computational cost. 
The best way to make a statistical sample of model galaxies 
is by combining large cosmological dark matt er simulations, 
such as the Millenni um Simulation (MS, Springel et al.l 
120051 ) or the MS-II (jBovlan-Kolchin et alJ [20091 ). with a 
method to include g a laxies . This approach was adopted by 
Ivan den Bosch et al. I (|2005l ). who used the conditional lu- 
minosity function technique that simultaneously optimizes 
the model match to the abundance and clustering of low 
luminosity galaxies in o rder to study the s atellite projected 
number density profile. IChen et al.l (|200fJ ) investigated the 
same statistic by assigning luminosities to dark matter 
structures so as to match the simulated cumulative circular 
velocity function to the SDSS cumulative galaxy LF. A 
si milar subhalo abund ance matching method was used 
by iBusha et al. I Hon]) to study the frequency of bright 
satellites around MW-like primaries. 

Semi-analytic models provide a more physically mo- 
tivated approach to including galaxies into dark matter 
simulations and have been show n to match a wi de va- 
riety of observational data (e.g iKauffmann et all 
Lacev et al.l 19931; Cole et al. 1994; IKauffmann et alj 



1993 



1999 



Somerville fe Primacklll999l; ICole et alj|200d ; [Benson et ail 
2002L iBaugh et alJ 120051 : iBower et al.l (2009 ; ICroton et ail 
200d : lGuo et al.ll201ll ). A large number of studies have used 
this techn ique applied to simu lations such as the MS-II, 
Aquarius (|Springel et alj |200ST ) and others to study vari- 
ous aspects of the galaxy population predicted in a ACDM 
model (e.glMunoz et aLl | 2009l; ICooper et al1l20ld; iLi et all 
20ld;lMacci6 et alJl20ld ; lFont et al.ll201ll ; IWang et alJl2oia 



Wang fe White! |2012|). In partic ular, the mock catalogues 



constructed by I Guo et all (l201lh were tested agai nst data 
from the SDSS in two studies. ISales et al.l (|2012r ) showed 
that the abundance of satellite galaxies as a function of pri- 
mary stellar mass in the SDSS DR7 spectroscopic catalogue 
was in good agr eement with this mode l. Considering fainter 
dwarf satellites, IWang fc White] |2012l ) studied the luminos- 
ity, colour distribution and stellar mass function using SDSS 
DR8 data, concluding that, apart from the model satellites 
becoming red too quickly when entering the halo of the pri- 
mary galaxy, many of the observed trends were reproduced 
in the model. 

In this paper, we will test the ACDM model 
and th e semi-analytic galaxy formation model GAL- 
FORM ijBower et alj|2006l) . by comparing the properties of 
model galactic satellite systems with those measured from 
the SDSS DR8 spectroscopic and photometric galaxy cata- 
logues. We will introduce a new series of GALFORM model 
galaxy catalogues based on MS-II that will allow us to ex- 
tend the predictions from the MS to less luminous galaxies. 

In Section 2 we describe the SDSS and model galaxy 
catalogue data that we use, and compare these two data 



sets. Section 3 contains summaries of the methods we use 
to select primaries and determine the satellite luminosity 
function and the projected number density profile. These 
two satellite galaxy distributions, upon which we will fo- 
cus in this paper, were calculated for the SDSS samples by 
Guo et al. (2011, hereafter Paper I) and Guo et al. (2012, 
hereafter Paper II). A verification that our estimators are 
unbiased is performed using the model galaxy samples. The 
results of the comparison between model satellite systems 
and those around similarly isolated SDSS primaries is pre- 
sented in Section 4. Implications for the model drawn from 
these comparisons are discussed in Section 5; we conclude 
in Section 6. Throughout the paper we assume a fiducial 
ACDM cosmological model with = 0.3, S1a = 0.7 and 
H = 70 km s _1 Mpc _1 . 



2 OBSERVED AND MODEL GALAXIES 

In this section we briefly review the data being used from 
the SDSS, which were described in more detail in Paper I 
and Paper II, before introducing the procedure to construct 
mock galaxy catalogues to compare with these data. 



2.1 SDSS galaxies 

Galaxies from both the spectroscopic and photometric sam- 
ples in the SDSS DR8 are used for this study. Isolated pri- 
mary galaxies, as defined in Paper I and Paper II, are se- 
lected from the spectroscopic survey, whereas satellites can 
come from either the spectroscopic or photometric surveys. 
With the relatively poorly constrained distances provided 
by the photometric redshifts, a statistical background sub- 
traction is performed to obtain an estimate of the satellite 
galaxy population around each of the primaries, as described 
in Section [3] below. 



2.2 Model galaxies 

The model galaxy catalogues were created using a combi- 
nation of large dark matter simulations to define the mass 
distribution and a semi-analytic model to place the galax- 
ies within this density field. E ither the Millennium Simula- 
tion (MS, ISpringel et alj|2005l) or the Millennium-II Simu- 
lation (MS-II. iBovlan-Kolchin et al.ll2009h was used to pro- 
vide the mass distribution. The former covers a large volume 
and hence contains many suitable isolated primaries, while 
the latter traces the mass distribution in a smaller volume, 
thus resolving structures containing lower luminosity galax- 
ies. While the two simulations trace the same number of 
particles (~ 10 10 ), the MS and MS-II simulation cubes are 
714 Mpc and 143 Mpc long respectively. The corresponding 



particle masses are ~ 1.23 x 10 y M and 9.9 x 10" M . 

The model galaxies populate the dark matter struc- 
tures according to the galaxy formation model GALFORM 
l|Bower et al.ll2006h . which includes reionization at high red- 
shift and energy injection from supernovae and stellar winds 
in order to prevent the overproduction of low luminosity 
galaxies. The luminosities of the most luminous galaxies are 
curbed by feedback from AGN. Parameters in the treatment 
of the galaxy formation processes have been chosen so as to 
produce as good a match as possible to the observed K band 
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Figure 1. The luminosity functions of all galaxies in the MS (blue 
dashed line) and MS-II (black solid line) cubes compared with 
observed luminosity functions of galaxies in the GAMA survey 
(the z < 0.1 subset from Loveday et al. 2012, red points) and 
SDSS (Blanton et al. 2005, green points). The observed luminosity 
functions are k-corrected (SDSS) or shifted (GAMA) to z = 0. 



LF of local galaxies. Fig. [T] shows the luminosity functions 
of all galaxies in the MS and MS-II simulation cubes. They 
match quite well with both the observed r band luminosity 
function of GAMA galaxies l|Lovedav et al.ll2012Ti shifted to 
z = by applying the offset r = 01 r — 0.22 (Blanto n et al.l 
120051 ). and the LF from the SDSS (jBlanton et al.ll2005t ). The 
drop in the MS LF at low luminosities reflects the resolu- 
tion limit, which corresponds roughly to M r = —16. Galax- 
ies placed into the MS-II should be complete significantly 
beyond this and thus suitable for studying faint satellites. 

For the MS, in addition to having the GAL- 
FORM model galaxies populating the simulation cube, flux- 
limited light-cone mock catalogues, either with or without 
the peculiar velocity included in the line-of-sight 'r edshift', 
have also been constructed (|Merson et alj |2012| ). These 
galaxies cover the redshift range z = 0.0 — 2.0. To simulate 
the photometric redshifts in the SDSS DR8 catalogue, we 
assign a photometric redshift with a random error to every 
galaxy that is fainter than m r = 17.7. 

Using both the MS and MS-II cubes of galaxies, we can 
test how robust our results are to the numerical resolution 
of the dark matter simulations. Comparison of the results 
obtained from the MS cube with those obtained from the 
light-cone mock catalogue, provides a test of the accuracy 
of our background removal and satellite distribution esti- 
mation procedures. Finally, the light-cone mock catalogue 
is intended to mimic the SDSS survey and provide a direct 
test of the model. 

When calculating scaled satellite number density pro- 
files, it is necessary to determine the value of r2oo (defined 
as the radius enclosing a mean total overdensity of 200 
times the critical cosmic value) associated with each pri- 
mary galaxy. Following Paper II, r2oo is estimated from the 
stellar mass, inferred from the luminosity and colour of the 
primary. This is converted to a hal o mass, M 2 00, us ing the 
abundance matching technique of iGuo et al.l (|2010l ). from 
which r2oo follows. The solid and dotted lines in Fig. [2]trace 
the distributions of r2oo estimated from the MS and SDSS 
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Figure 2. The distribution of estimated T200 values for model 
primary galaxies (solid lines), real r^oo values for their associ- 
ated parent haloes (dashed lines) and estimated T2oo values for 
primaries selected from the SDSS (dotted lines). From top to bot- 
tom, the panels show results for primary magnitude bins centred 
on M c = —21.0, —22.0 and —23.0 respectively. The vertical lines 
are the means of the corresponding distributions. 



primaries respectively, with the different panels showing re- 
sults for different luminosity primaries. Vertical lines show 
the mean values of the distributions, which differ by no more 
than about 15 per cent in all cases. This similarity between 
estimated satellite system sizes in the mock and SDSS sug- 
gests that scaling the satellite number density profiles by 



© 0000 RAS, MNRAS 000, 000-000 



4 Guo et al 





Figure 3. Colour-magnitude diagrams of spectroscopic galaxies in SDSS (top left) and of galaxies in the MS mock light-cone (top right) 
and the two different simulation cubes (MS, bottom left and MS-II, bottom right). Contours represent lines of constant galaxy number 
density and the straight lines indicate the colour cuts used to define red and blue primaries. The same cut is used for both the MS and 
MS-II. 



the system size should not create any large systematic dif- 
ferences between the real and mock results. 

The distributions of real r2oo values, inferred from the 
dark matter distribution in the simulation cube, are shown 
with dashed lines in Fig. [2] For the lower luminosity pri- 
maries, the real and estimated V2oo distributions are sim- 
ilar. However, in the bottom panel of Fig. [5] it is appar- 
ent that there is a significant population of physically small 
haloes surrounding primary galaxies whose stellar mass cor- 
responds to a larger halo. As pointed out in Paper II, this 
is probably due to the fact that the stellar mass varies little 
with halo mass for these most luminous systems. Thus, any 
small scatter in stellar mass gives rise to a large change in the 
value of r2oo inferred from abundance matching. This is an 
important source of potential bias when trying to measure 
concentrations of satellite distributions around the most lu- 
minous primaries. In this case, a concentration derived using 
the value of r2oo estimated from stellar masses will not nec- 
essarily reflect the concentration as defined with respect to 
the halo mass. 

Finally, as we will be investigating the colour depen- 
dence of the results, in Fig. [3] we compare the colour dis- 
tributions of galaxies in the SDSS and model galaxy cat- 
alogues. The distribution of SDSS galaxies in the colour- 



magnitude diagram is shifted along the g — r axis com- 
pared to that of the model galaxies. Thus, while we choose 
°'°(S - r)cS SS = °- 15 ~ 0.024M r as the line dividing the red 
and blue populations in the SDSS this boundary is placed 
at (3-r)^t' MS " n = -0.28-0.04A'/,. for the model galaxy 
populations. These cuts are shown in the 4 panels of Fig. [3] 
where the effect on the colour magnitude distribution of in- 
cluding the survey flux limit can be seen for the MS. This 
removes the bulk of the low luminosity galaxies that are 
present in the z = snapshot and makes the resulting light- 
cone galaxy population significantly more like that in the 
SDSS. The fraction of the -21 < M r < -19 model light- 
cone galaxies that are defined as red is 0.51, the same as 
that for the SDSS. Had we adopted the SDSS colour cut, 
the model would only have had a red galaxy fraction of 0.28. 
For galaxies in the magnitude range —20 < M r < —19, the 
red fractions are 0.46 and 0.38 for the model and SDSS re- 
spectively, when using the cuts shown in Fig. [3] Thus, apart 
from a magnitude-dependent shift in the colours of galaxies, 
which we can correct for with the different colour cuts, the 
global red galaxy fractions are quite similar between model 
and SDSS data sets at magnitudes that we will be consider- 
ing for the satellite galaxies. 

The MS and MS-II simulations, once populated with 
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galaxies according to the semi-analytic model, contain very 
similar galaxy populations. Thus, it is appropriate to use 
them interchangeably depending upon which is more impor- 
tant: having a large volume containing many primaries, or 
being able to resolve low luminosity satellite galaxies. How- 
ever, if we want to compare results from the MS-II cube 
of galaxies with those from the SDSS flux-limited survey, 
then we still need to demonstrate, using the MS, that our 
methods recover from the light-cone mock surveys the same 
satellite distributions as are present in the simulation cubes. 



3 METHOD 

In this section we briefly review the procedure used to deter- 
mine the satellite LF and projected number density profile 
for SDSS, described more fully in Papers I and II, before 
detailing how these distributions are determined from the 
various different types of model galaxy catalogue. The qual- 
ity of the recovery of these distributions is quantified by 
comparing those determined from the MS cube of galaxies 
with those from the flux-limited mock light-cone surveys. 

Primary galaxies are selected to have spectroscopic 
redshifts in the SDSS and to have magnitudes, M p , sat- 
isfying M c - AM bin < M p < M c + AM bin . We choose 
M c = -21.0, -22.0, -23.0 and AM bin = 0.5. Further, the 
primaries should be isolated in the sense that no other 
galaxy within AMf a j n t = 0.5 magnitudes lies within a pro- 
jected distance of 2i?i nncr and is sufficiently close in red- 
shift. 'Sufficiently close' is defined as a difference in spec- 
troscopic redshift of less than Az s = 0.002 or, for galaxies 
without a spectroscopic redshift, with a photometric red- 
shift within aptjp, where ap — 2.5 and op is the pho- 
tometric redshift error defined in Paper II. R mn er repre- 
sents the projected radius within which satellites may re- 
side, and the variable -Router defines the outer edge of an 
annulus within which the local background is estimated. 
We adopt the same values of (Rinner, Router) as in Paper II: 
(0.3, 0.6), (0.4, 0.8) and (0.55, 0.9) Mpc for primaries in mag- 
nitude bins M c = -21.0, -22.0 and -23.0 respectively. Only 
galaxies brighter than m r lm = 20.5 are considered. 

Only sufficiently close galaxies in redshift are included 
when counting the potential satellites within the pro- 
jected radius Rinner and making the local background es- 
timate from the surrounding annulus out to Router- The 
background-subtracted satellite systems are stacked for pri- 
maries in each absolute magnitude bin to provide estimates 
for the mean satellite LFs and projected number density pro- 
files of satellites more luminous than a particular absolute 
magnitude, as described more fully in Papers I & II. 

The procedure described above is applied to the SDSS 
itself and also to the MS redshift space light-cone mock cat- 
alogue. However, different estimation procedures are used 
for the light-cone mock with real space (rather than redshift 
space) galaxy coordinates and the galaxy populations in the 
simulation cubes. With the real space positions, it is possible 
to define isolated primaries as having no bright neighbours 
within a sphere of radius 2Ri nri cr- The satellites within a 
sphere of radius Rmner can also be determined without any 
need for background subtraction. Using only the real space 
satellites, the estimation of the satellite luminosity function 



is straightforward. For the jth magnitude bin, the average 
satellite LF is estimated using 

-^rreal sat _ ^jyreal sat ^jy^ 

J ^yreal prim ' * ' 

where A r 1 rcal sat is the number of satellites around primary 
i and A^ eal prlm is the number of primaries contributing to 
the jth bin of the LF. 

While we need no correction for interlopers in real space, 
the process of estimating the mean projected number den- 
sity of satellite galaxies is such that a fair comparison with 
light-cone data still requires the subtraction of a background 
estimated from an outer area. Hence, the projected number 
density profiles of satellites brighter than M tlun in real space 
are determined using all galaxies within a cylinder of pro- 
jected radius Router and length 2R ou tcr, centred on the pri- 
mary. These galaxies are projected onto a plane and provide 
the sources for the potential satellites/background for pro- 
jected radii less/greater than Rmner. The projected number 
density profile is determined using 

where iVy is the number of galaxies brighter than M run 
within a projected distance Rinner of the ith primary and 
in the jth projected annulus, and N^ ck is the corresponding 
number of galaxies in the projected outer annulus, Rinner < 
r < Router- A\. is the area contributed by the ith primary to 

the jth annulus for the detection of satellites brighter than 

M tnm 

{A irtrun - J/flim 

M trun > M] im , (i) 

where Aij is the area of the jth annulus surrounding the 
ith primary and A/' lm is the absolute magnitude that corre- 
sponds to the apparent magnitude limit of the mock cata- 
logue. A° utCT is the corresponding area in the outer annulus 
surrounding the ith primary. 

The comparison of the projected number density pro- 
files, before and after background subtraction, with that 
formed by simply projecting the galaxies within a sphere of 
radius Rmner of the primary galaxy is shown in Fig. [4] The 
results indicate that the projected profile after subtracting 
the background very accurately recovers that estimated di- 
rectly from the galaxies within the inner area (r < Rmner). 
The impact of the background subtraction is small and lim- 
ited to radii near to Rinner- This establishes that the method 
for calculating the background subtracted projected number 
density profile from a real space light-cone survey provides 
an unbiased estimate of that produced when only satellites 
within a 3D distance Rinner are used. We can now compare 
these real space profiles with those from redshift space light- 
cones. 

Having described the three different types of model 
galaxy catalogues made using the MS and how the satel- 
lite distributions are determined from each of them, we are 
now in a position to test the accuracy of our estimation pro- 
cedure. Fig. [5] shows the satellite LF and projected number 
density profiles estimated from each of the model galaxy cat- 
alogues: simulation cube, real space light-cone and redshift 
space light-cone. The agreement between the two different 
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Figure 4. The comparison of real space projected number den- 
sity profiles estimated by different methods. The black line is the 
profile estimated by projecting the satellites within spheres of ra- 
dius -Ri nn er around primaries. The dashed blue curve shows the 
projected number density profile of all galaxies within a line-of- 
sight distance i? uter of a primary. The solid blue curve is the 
background-subtracted case, where the background is estimated 
from the outer annulus with -f?i nncr < r < -Router • Vertical dashed 
lines show Ri nnC r and Router. 



satellite distributions measured from all three catalogues at 
Af c = —21.0, —22.0 and —23.0 provides strong support that 
our results are unbiased, and that our technique for back- 
ground subtraction is appropriate. Given that the isolated 
primaries are a specially selected subset of the differentially 
clustered galaxy population, one would expect that a local, 
rather than global, background subtraction would be ap- 
propriate. The fact that the results from the MS light-cone 
mock match well with those from the MS cube of galaxies 
suggests that we can use results from the MS-II simulation 
cube at z — when comparing with low luminosity galaxies 
in the SDSS DR8. 
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Figure 5. Satellite LFs (top) and number density profiles (bot- 
tom) estimated from light-cone mock catalogues and directly from 
the MS cube. The results for primary magnitude bins centred at 
M c = —21.0, —22.0, —23.0 are shown in black, blue and red re- 
spectively. Solid and dashed lines correspond to results from red- 
shift and real space light-cones respectively, whereas the points 
show the results for the whole volume in the simulation cube. 



4 RESULTS 

Having verified that the model galaxies have similar distri- 
butions of luminosity and colour to those in the SDSS and 
that our local background subtraction procedure produces 
unbiased estimates of the satellite LF and projected num- 
ber density profile, we now compare the model and observed 
satellite systems in more detail. 

4.1 Dependence on primary and satellite 
luminosity 

The top panel of Fig. [6] shows the satellite LFs for the pri- 
mary magnitude bins M c — —21.0, —22.0, —23.0 estimated 
from both the SDSS and MS-II mock data. For AM r < 5, 
the model and SDSS satellite luminosity functions generally 
agree well. However, there is a steepening of the LF for lower 
luminosity satellite galaxies in the SDSS that is not present 
in the model. The MS-II was used for this comparison, so 
a lack of numerical resolution should not be responsible for 
this deficit. 

The middle panel of Fig. [6] shows the projected number 



density profiles of satellites brighter than Ml™ = —19.0 for 
the different primary samples. This limit is chosen, following 
Paper II, to ensure that the profiles for different luminosity 
bins are measured from a large enough sample. In all cases 
the profile approximately follows E(r) oc r -1 ' 5 , with the 
amplitude reflecting the fact that more luminous primaries 
host more satellites. For most radii, the model reproduces 
the amplitude observed in SDSS. In detail though, there is 
a factor 2 difference at ~ 30 kpc for both M c — —21.0 and 
—23.0. The model underproduces satellites at this radius for 
the least luminous primaries and overproduces them for the 
most luminous primaries. 

Once rescaled in both radius and projected number den- 
sity, as shown in the bottom panel of Fig. [6] the model 
and SDSS profiles line up well outside the region at 0.05 < 
r/r2oo <, 0.1. As noted in Paper II, the slight deficit of satel- 
lite galaxies in the inner regions of the SDSS for Af c = —23.0 
primaries, relative to the model, may reflect difficulties iden- 
tifying low luminosity satellites in regi ons where the back - 
ground light subtraction is significant l)Aihara et al.ll20Tll ). 
This had the effect of slightly decreasing the fitted concen- 
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AM r log(r/Mpc) 





Figure 6. Model (MS-II, solid lines) and SDSS (points) satel- 
lite LFs (top), projected number density profiles (middle) and 
normalized profiles (bottom). The results for primary luminosity 
bins M c = —21.0, —22.0, —23.0 are shown in black, blue and red, 
respectively. 



tration of satellites around the most luminous primaries rel- 
ative to the low luminosity primary bins. 

The dependence of the satellite projected number den- 
sity profile on satellite luminosity for primaries with —22.5 < 
M p < —21.5 is shown in Fig. [7] Satellites are split into two 
bands of luminosity: —21.0 < M s < —19.0, representing 
the objects that contributed to the profile in Fig. H3 and 
— 19.0 < M s < —17.0, showing the behaviour of lower lumi- 
nosity satellites. The top and middle panels show the satel- 



Figure 7. The dependence of the scaled and unsealed satellite 
density profiles on satellite luminosity for primaries in the range 
—22.5 < M p < —21.5. Unsealed (top) and scaled (middle) profiles 
of satellites are shown for two different luminosity bands: —21.0 < 
M s < -19.0 (black) and -19.0 < M s < -17.0 (blue). The solid 
lines and points correspond to results from the model and SDSS 
respectively. The bottom panel shows scaled profiles for lower 
luminosity satellites in the model. 

lite projected number density profiles before and after scal- 
ing respectively. While generally in good agreement, there is 
an excess of lower luminosity satellites in the SDSS relative 
to the model at ~ 30 kpc. 

The 'bump' in the projected number density profile in 
the SDSS data is present only for the lower luminosity satel- 
lites. One is therefore tempted to ask if it is present at any 
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satellite luminosity in the mock catalogues. This question is 
answered in the bottom panel of Fig. [7J where no compara- 
ble deviations are seen for satellites with M s < —13, which 
have indistinguishable profiles from those of brighter satel- 
lites. The profile for satellites with — 13 < M s < — 11 shows 
a slight change in shape over a large range of scales relative 
to the other sets of satellites, but nothing as pronounced as 
is seen in the SDSS. 



4.2 Dependence on primary and satellite colour 

As a further test of the galaxy formation model, which has 
so far largely succeeded in reproducing the satellite LF and 
projected number density distributions, we now split the 
primaries and satellites by colour using the slightly different 
colour cuts described in Section ^. 2l for the model and SDSS. 

Fig. [8] shows how the satellite luminosity function and 
projected number density profile depend upon the colour 
of the primary around which they are being measured, and 
compares the results from the model and the SDSS for pri- 
maries with —22.5 < M p < —21.5. For the adopted colour 
cuts, the vast majority of primaries are classified as red in 
both the model and SDSS. Given that when not split by 
colour the model produces a good match to these satellite 
distributions, it is no surprise that the satellite populations 
around red primaries agree well between model and obser- 
vations. However, the blue primaries in the mock catalogues 
are deficient in satellites by a factor of 2 — 3. The excess satel- 
lites around SDSS primaries span the range of luminosities 
and radii being considered here, with a slight tendency to be 
at lower luminosity and nearer to the primary than for the 
satellites present around the model primaries. The scaled 
profiles in the bottom panel of Fig. [8] show that satellites 
around blue SDSS primaries are slightly more concentrated 
than those around model primaries. 

In the top two panels of Fig. [8] results for both the MS 
and MS-II are shown. Because of its limited resolution the 
LF in the MS becomes incomplete at AM ~ 5.5, but the 
LF in the MS-II is well resolved down to AM ~ 7. On the 
other hand, in the smaller volume of the MS-II, there is a 
relatively small number of primaries and, as a result, the 
projected number density profile of satellites brighter than 
M a = — 19 is noisy (and it is therefore omitted in the lower 
panel of the figure). In the regions where both simulations 
are well resolved and sampled, their results are consistent. 

The satellites themselves can be divided into red and 
blue subsets and their distributions around primaries with 
—22.5 < M p < —21.5 are shown in Fig. [5] Once again both 
MS and MS-II results are shown in the top two panels, while 
the noisy MS-II results are not presented for the scaled pro- 
jected number density profiles in the bottom panel. The fact 
that the colour cuts in the model and SDSS samples yield 
completely different red and blue fractions is immediately 
apparent in the satellite LFs, with the model satellite sys- 
tems dominated by red satellites to a coincident ally simi- 
lar extent as the blue satellites dominate around SDSS pri- 
maries. For the most luminous satellites, the SDSS blue and 
red fractions converge, whereas this does not happen for the 
model. 

The shape of the projected density profiles of red satel- 
lites are very similar for the model and SDSS systems, once 
the difference in amplitude has been scaled away, as shown in 
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Figure 8. The dependence on primary galaxy colour of the satel- 
lite LF (top), unsealed (middle) and scaled (bottom) number den- 
sity profiles for primary galaxies of magnitude —22.5 < M p < 
—21.5. Solid and dashed lines show model results for primaries in 
the MS and MS-II respectively, whereas the points are for SDSS. 
All profiles are for satellites more luminous than M B = —19.0. 



the bottom panel of Fig. [9] However, this is not the case for 
the blue satellites, which have a much less concentrated dis- 
tribution around model primaries. For real systems of blue 
satellites, the distribution is only slightly less concentrated 
than it is for the red satellites. 
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Figure 9. Dependence of satellite distributions on the colour 
of the satellite galaxy. Symbols and line types have the same 
meaning as in Fig. [8] 



5 DISCUSSION 

The comparison of galactic satellite systems in the model 
with those in the SDSS shows that the dependence of the 
satellite distributions on primary and satellite luminosity 
is captured quite well by the GALFORM model. This is a 
non-trivial success of the model since its parameters were 
adjusted merely to match the global K-band LF of galaxies, 
with no direct reference to satellite systems. There is an ex- 
cess of very low luminosity satellites around SDSS primaries 
relative to the model, and the projected number density pro- 



files are up to a factor 2 discrepant within ~ 30 kpc, but the 
agre ement is general ly good. 

iTal et al.l (|2012n also studied the radial distribution of 
satellite systems around bright primary galaxies using SDSS 
data. Their primaries were LRGs at 0.28 < z < 0.40 with no 
isolation criteria applied and hence often in groups of bright 
galaxies, and thus are different to those studied here in a 
few respects that may well be important. They found the 
projected number densit y was well fitted by a comb ination of 
a projected NFW profile (jNavarro et al.lll99fj| . ll997h for large 
radii and a central steeper profile that follows the stacked 
light profile of the LRGs. This central bump is similar to 
that seen for the lower luminosity satellites around primary 
galaxies shown in Fig. [7J but is, in contrast, most apparent 
in the high luminosity satellites around the LRGs. 

The differences between the model and SDSS satel- 
lite systems are greater when the colour of the satellites is 
considered. Even the distribution of galaxies in the colour- 
magnitude plane shows that the model has too high a frac- 
tion of low luminosity red galaxies relative to the SDSS. The 
model blue satellites are both significantly depleted and very 
much less concentrated relative to either blue or red SDSS 
satellites, which have a projected number density profile like 
that of the red model satellites. These pieces of evidence 
point to the model being too ready to convert low luminos- 
ity blue galaxies to red ones . 

IWeinmann et all (|2012f ) suggest that galaxy formation 
models have a generic problem with maintaining sufficient 
gas to form stars in lower mass galaxies at low redshift. The 
satellite galaxies we ar e considering are somew hat larger 
than those studied by IWeinmann et all ((2012) , and our 
choice of different colour cuts for defining red and blue galax- 
ies in the model and SDSS samples reduces global systematic 
colour differences. For instance, the overall fraction of blue 
galaxies in the model, at the magnitudes of the satellites 
that we focus on, is very similar to that in the SDSS. Satel- 
lite galaxies in our study constitute only a small fraction of 
the total population of galaxies. Thus, a dramatic change in 
the satellite properties will leave very little imprint on the 
global LF. 

Alternatively, it could be that semi-analytical galaxies 
turn red too rapidly after accretion into larg er haloes. This 
idea was investigated by iFont et all |2008), who changed 
the GALFORM treatment of gas stripping from subhaloes 
as they enter the virial radius of large primary galaxies. 
Rather than hot gas being stripped from subhaloes imme- 
diately as they enter th e virial radi u s, a m ore gradual loss 
of gas is adopted in the IFont et al.l i|200Sft model. This al- 
lows a relatively extended period of star formation to occur 
and the possiblity of bluer satellites. We have performed 
our analysis on a model galaxy population construc ted us- 
ing this particular variant of the lBower et al.l (120061 ) GAL- 
FORM model. While the typical colours of galaxies do be- 
come slight ly bluer, the num ber of blue satellites per pri- 
mary m the lFont et al. (2008) model increases only slightly, 
as can be seen in Fig. 1101 The shape of the blue satellite 
profile improves significantly, with the extra blue satellites 
being found preferentially at small radii. However, the abun- 
dance of red satellites is also increased by this modification 
to the GALFORM model, because satellites generally be- 
come more luminous as a result of the more extended pe - 
riod of star formation. As a result, the IFont et all (|2008T ) 
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Figure 10. The effect on the model satellite LF (top) and num- 
ber density profiles (lower two panels) of changing the treatment 
of hot gas stripping in GALFORM for primary galaxies of mag- 
nitude —22.5 < M p < —21.5. The points are the results mea- 
sured from SDSS data, whereas solid and dashed lines are the 
results measured from catalogues constructed using the Bower et 
al. (2006) and Font et al. (2008) GALFORM models respectively. 
The middle panel shows the profiles of all satellites, whereas the 
lower panel shows them split by colour. 



model overproduces satellite abundances overall, as shown 
in the upper two panels of Fig. 1101 The overproduction is 
most discrepant with the data at low luminosities. 

Many important astrophysical processes combine to de- 
termine the distribution of low luminosity galaxies. There- 
fore the distributions of satellite galaxies will depend sensi- 



tively on aspects of the galaxy formation model. Given that 
the treatment of gas stripping can have the large impact 
shown in Fig. 1101 o ne is dr a wn to conclude that the ability 
of the default Bow er et all (2006) model to match the to- 
tal satellite LF and projected number density profile of the 
SDSS systems was far from inevitable. 



6 CONCLUSIONS 

Using model galaxy catalogues constructed using large dark 
matter simulations and a semi-analytic galaxy formation 
model, we have tested the accuracy of our procedures for 
measuring properties of the distribution of satellite galaxies 
around bright, isolated primary galaxies. We find that our 
local estimation of the abundance of background galaxies 
yields unbiased estimates of the satellite galaxy luminosity 
function and projected number density profile. The agree- 
ment between results in the MS and MS-II galaxy catalogues 
in their region of overlap shows that our results are numer- 
ically converged and allows us to extend the dynamic range 
of the model predictions. 

Comparing the model predictions with those measured 
for satellite systems in the SDSS, we find that the de- 
pendence of the satellite LF is matched well for M c = 
-21.0, -22.0, -23.0 and AM r < 5. Lower luminosity satel- 
lites are increasingly underpredicted by the model. The pro- 
jected number density profile is also well reproduced at radii 
greater than ~ 30 kpc. At smaller radii, deviations in the 
abundance by a factor of two are apparent. These differ- 
ences between model and SDSS are seen most strongly in 
the low luminosity satellites, which show an excess in the 
SDSS relative to the extrapolation of the power law from 
larger radii, which describes the inner regions of the model 
satellite systems. 

Splitting the sample into red and blue galaxies pro- 
duces more dramatic differences between the model and 
SDSS results. The model places a factor 2 — 3 fewer satellites 
around blue primaries than are present around comparable 
SDSS primaries. However, the discrepancy between model 
and SDSS is even larger when considering the colours of 
the satellites. The model satellites are predominantly red, in 
contrast to the blue-dominated SDSS satellite galaxy pop- 
ulation. Furthermore, what model blue satellites there are 
have a significantly more extended distribution around their 
primary galaxy than is seen for either the SDSS blue satel- 
lites or the red satellites in the model and SDSS. 

The generally successful comparison of the GAL- 
FORM model with the SDSS data performed here provides 
a non-trivial validation of the assumptions and framework 
of this kind of modelling. At the same time, the failure of 
the model to account for the observed colour dependence of 
the satellite properties demonstrates that the model is in- 
complete and that important physical processes, almost cer- 
tainly related to the rapidity with which infalling satellite 
galaxies turn red, are not being faithfully modelled. Since 
a similar sho r tcomi ng is present in the independent model 
of iGuo et al.l (|201ll ). this problem seems deep-rooted and is 
worthy of further investigation. 
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